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REDUCED MODELS FOR LINEARLY ELASTIC THIN FILMS 
ALLOWING FOR FRACTURE, DEBONDING OR DELAMINATION 

JEAN-FRANgOIS BABADJIAN, DUVAN HENAO 


Abstract. This work is devoted so show the appearance of different cracking modes 
in linearly elastic thin film systems by means of an asymptotic analysis as the thickness 
tends to zero. By superposing two thin plates, and upon suitable scaling law assumptions 
on the elasticity and fracture parameters, it is proven that either debonding or transverse 
cracks can emerge in the limit. A model coupling debonding, transverse cracks and 
delamination is also discussed. 
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1. Introduction 

It is experimentaly observed that thin films systems can essentially develop two different 
crack patterns: either transverse cracks channeling through the thickness of the film, or 
planar debonding at the interface of two layers. In classical fracture mechanics, a threshold 
criterion on the energy release rate drives the propagation of a crack along a prescribed 
path. Within this framework, m described different possibilities of failure modes. In 
|42] . a reduced two-dimensional model of a thin film system on an elastic foundation is 
proposed, and the propagation of different crack modes is discussed. This model is later 

Key words and phrases. Fracture mechanics. Functions of bounded deformation, F-convergence. 
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recast as an energy minimization problem, based on the variational approach to fracture 
of [13j . hrst in [32] in a simplihed one-dimensional setting where transverse cracks are 
represented by a hnite number of discontinuity points for the displacement, then in [33| in 
the full two-dimensional case. The total energy is the sum of a bulk energy (including the 
elastic energy of the him outside the transverse cracks) and a surface energy of Griffith type 
(including the area of the transverse fractures and of the debonded regions). Concretely, 
it takes the form 


Eiu,T,A)= [ 

J Uj\ 


\r L-^/ + 2h/ 


eaa{u)epj3{u) + fj.fea/3{u)eap{u) 


dx 


Ae(u):e{u) 



u — w\‘^ dx + Kf length(r) -|- Kf, area(A); 


we proceed to explain each term separately. The region a; C denotes the basis of a thin 
him = uj X {0,hf) bonded on a inhnitely rigid substrate, where hf is the thickness of 

the him and e = L = diamcj, is a non-dimensional small parameter. The transverse 
cracks are of the form T x (0,/ij) where T is a one-dimensional object which can be 
thought of as the union of a hnite number of closed curves, which are themselves part 
of the unknowns of the problem. The delamination zone A C a; is also an unknown. 
The in-plane displacement of the him at the interface with the substrate is denoted by 
u : UJ \r —?> M^. The fracture toughness kj is a material property of the him, while Hb 
measures the strength of the bonding between the him and the substrate. The reduced 
linearly elastic energy Ae{u) : e{u) is well-known and rigorously derived in the Kirchhoff- 
Love theory of elastic plates m- 

It remains to explain the term := ^ fu;\A 1'“ ~ w\‘^dx. The map m : u; —>■ is 

given and represents the displacement of the substrate. Since the substrate is assumed 
to be inhnitely rigid, w is the same displacement it would undergo if the him were not 
present. Note also that we are only considering planar displacements of the substrate. 
The energy Ec represents the price to pay in order for the him to deform differently 
from the substrate. It only has to be paid in cj \ A because in A the him is no longer 
attached to the substrate. By regarding the him and the substrate as a single elastic 
body, it is seen that Ec has the form dx of a Barenblatt’s cohesive-zone surface 

energy [9], where [u] represents the jump of the displacement across the debonding zone, 
in this case with the integrand g{s) = These cohesive energies are considered, in 

particular, in the existing analytical studies of delamination problems, e.g. [niETlETlEHI. 
Apart from fracture mechanics, this type of integrals also appear in the study of Winkler 
foundations m, with applications as varied as the understanding of the seismic response 
of piers, chromosome function, or the mechanical response of carbon nanotubes embedded 
in elastic media (see [34| and the references therein). In this setting, Ec is interpreted as 
the effective energy of an elastic foundation, understood as a continuous bed of mutually 
independent, linear, elastic springs, hence the appearance of a reaction force of the form 
gb{u(x) — w{x)) (corresponding to the quadratic energy ^\u — top) in response to the 
relative displacement of the body supported on the foundation. 

The question addressed by this paper is the rigorous derivation of E{u, T, A) from 
three-dimensional linearized elasticity in the limit as e —)• 0. This involves 
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(1) the derivation of a reduced Griffith model for the initiation and propagation of 
cracks in a thin film, and, 

(2) the justification of the cohesive energy for the debonding at the interface. 

Previous studies of the first problem include [IH [121 [H |33] , which consider scalar-valued 
problems or generalizations which are incompatible with 3D linear elasticity (coercivity 
assumptions of the form W{F) > C{\FY‘ — 1) are imposed on the stored-energy densities 
VP(Vu)), and [26], where the thin film is linearly elastic but the path and the geometry 
of the crack are specified a priori (it has to be of the form P x { 0 ,hf) and it can only be a 
single crack). In Theorem 15.II below, we present the first complete result for the reduction 
of dimension of a brittle linearly elastic thin film, without any prior assumption on the 
geometry or the topology of the cracks. As is well-known m, this falls in the framework 
of free discontinuity problems, where a satisfactory mathematical treatment can be done 
in the space of special functions of bounded deformation. We adopt usual scalings for 
the elastic and fracture parameters, and show the convergence to a reduced model where 
admissible cracks are vertical, and admissible displacements have a Kirchhoff-Love type 
structure (the out-of-plane displacement is planar, while the in-plane displacement is affine 
with respect to the out-of-plane variable). The main difficulty is to establish a compactness 
result on minimizing sequences (Propositions 15.11 and 15.2p showing the structure of limit 
displacements and cracks with finite energy. It uses tools of geometric measure theory and 
fine properties of bounded deformation functions. 

The justification of the cohesive energy Ec is also very delicate. In [29l [211 [211 [H] a 
free discontinuity model with a cohesive fracture energy is obtained as the P-limit of an 
Ambrosio-Tortorelli functional in which the constraint z > ^/£ is imposed on the internal 
damage variable, £ being the width of the damage zones. This is in the spirit of considering 
the possibility that what macroscopically would be regarded as fracture is actually a strain 
mismatch that is continuously accommodated through a very thin layer of a very compliant 
material. However, it is unclear whether their approach is suitable for the study of thin 
films, in particular if it could explain that only transverse fracture and planar debonding at 
the interface can be observed in the limit. Cohesive-type energies have also been obtained 
by homogenization in mm as the limit of a Neumann sieve, debonding being regarded as 
the effect of the interaction of two films through a suitably periodically distributed contact 
zone. A different derivation of a cohesive fracture energy (albeit with a positive activation 
threshold) can also be found in [T5|, in this case as a result of the homogenization of brittle 
composites with soft inclusions. 

Here we consider the problem of deriving Ec based on conceiving the interface between 
a bimaterial system as a very thin layer of a third phase, occupying the region Hfe = 
00 X [— /if,,0], where hf, represents its thickness. The expectation is to recover the cohesive 
delamination energy from the elastic energy of the bonding layer when it is made of a 
material that is increasingly more compliant as e —)• 0. A scaling law is then proposed for 
the Lame moduli of the adhesive, of the form (A^,//^) = £'^{Xh,Hb)i for some fixed Xb, Hb 
and some exponent q to be determined. Calling £b '■= ^ to the aspect ratio in Xlb and 
using the rescaled displacements Ua{x\x^) := Va{Lx', hbX^), U 3 (x',X 3 ) = hbV^lLx', X 3 ), 
where x' G oj/L, x^ G [—1,0] are non-dimensional rescaled spatial variables, we are able 
to write the energy of the bonding layer in the form 


4 


JEAN-FRANgOIS BABADJIAN, DUVAN HENAO 


with 


Je{v) ■= \ ! {4 

eaa{u)eijp{u) + 2^ibeaii{u)eap{u) 




+ 


2Xbeaa{u)e33{u) + 4:fJ.bea3{u)ea3{u) 


1 



(Afe + 2fib)e33{u)e33{u 


)} 


dx. 


If Je remains bounded as e — > 0, due to the in front of the third term, U 3 is expected 
to be planar in the limit; if the displacement of the substrate is planar, this means U 3 = 0 
outside the delamination zone. On the other hand, due to the coefficient for the 
first term, we expect the in-plane gradient to be irrelevant. Thus, the bonding layer 
is expected to behave according to ^ 9 aU 3 daU 3 , which is minimized if the strain 
mismatch between the film and the substrate is accomodated by an affine transition in the 
X 3 variable, giving rise to the cohesive energy Ec- The assumption that Jg is bounded in 
the asymptotic analysis corresponds to the energy in the bonding layer being of the same 
order of magnitude as the elastic energy of the film, which scales as hf. This yields the 
scaling 

Without loss of generality, in this paper we consider the case when the thicknesses of 
the film and of the bonding layer have the same order of magnitude and q = 2 , that is. 
For the effect of other scaling assumptions, we refer to [M] . 

The above heuristics were made rigorous in [33] in the simplihed case of scalar displace¬ 
ments. In the anti-plane case, where the problem becomes scalar, a simple adaptation of 
that result shows the convergence to a model coupling transverse cracks, cohesive transi¬ 
tions as long as the in-plane displacement is below a precise threshold, and delamination 
when the threshold is overpassed (Theorem [6T|) . In the full vectorial linearly elastic case, 
the reduced model was rigorously derived in |34j for the problem of Winkler foundations, 
that is, when the displacements are Sobolev maps so that neither the film nor the bonding 
layer are allowed to undergo fracture. In Theorem 14.11 below we give a simpler proof of 
the same result. 

We have been unable to prove the convergence to E{u,T, A) in the case of interest 
of a linearly elastic bonding layer which may undergo fracture. We limit ourselves to 
present some partial results which, in our opinion, ought to be considered in any attempt 
to establish the desired T-convergence. We prove an energy upper bound by constructing, 
for every admissible limit displacement, an optimal recovery sequence (Proposition 16. ij) . 
What remains open is to establish the optimality of the affine transitions in the X 3 variable 
in order to accommodate the mismatch between the film and the substrate. Indeed, the 
ability to break gives the bonding layer the opportunity to reduce its elastic energy by 
performing a periodic sequence of small rotations (Example 16.3p . This implies that the 
delamination zone cannot be identihed just by taking the orthogonal projections of the 
jump set of the displacement, as is done in the Sobolev and scalar cases. As a possible 
remedy, we consider instead, “almost vertical” projections. We are able to prove a surface 
energy lower bound (although with a bad multiplicative constant) and to show the validity 
of the desired bulk energy lower bound under the assumption that the minimizing sequence 
satishes better a priori estimates than just the energy bound (Lemma 16. 7p . 
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We end this Introduction by mentioning [36] , where a Griffith energy for the debonding 
at the interface is obtained as the limit elastic energy of a thin bonding layer in a problem 
involving a damage internal variable. The techniques of that paper may prove relevant in 
the derivation of the reduced model E(u,T, A) discussed in this paper. 

The paper is organized as follows: Section [2] is devoted to introduce various notations 
used throughout this work. In Section [S] we precisely decribe the model and perform a 
scaling to make the problem more tractable from a mathematical point of view. Section 
m investigates the asymptotic analysis in the absence of cracks, and evidences the appear¬ 
ance of a debonding type limiting energy fTheorem 14.ip . In Section [S] we carry out the 
analysis a linearly elastic thin film, and show the emergence of transverse cracks (Theorem 
EH). Finally, Section [6| discusses the interplay between transverse cracks, debonding, and 
delamination. 


2. Notation and preliminaries 

If a and b € M”, we write a-b = Euclidean scalar product, and we denote 

the norm by |a| = y/a ■ a. The open ball of center x and radius g is denoted by Bg{x). If 
X = 0, we simply write Bg instead of Bg{0). 

We denote by the set of real m x n matrices, and by the set of all real 

symmetric n X n matrices. Given two matrices A and B £ we let A : i? := iv{A^B) 

for the Frobenius scalar product, and |^| := for the associated norm {A^ is 

the transpose of A, and tr(A) is its trace). We recall that for any two vectors a £ and 
b £ ffi”, 0 ( 8)6 £ stands for the tensor product, i.e., {a®b)ij = Uibj for all 1 < i < m 

and 1 < j < n. If m = n, then aQ b := ^(a( 8)6 + 6 ( 8 )a) £ Mgy^ denotes the symmetric 
tensor product. 

Given an open subset U of M” and a hnite dimensional Euclidean space X. We 
use standard notations for Lebesgues spaces U’{U;X) and Sobolev spaces H^{U\X) or 
IT^’P(C/; X). We denote by Xi{U;X) the space of all X-valued Radon measures with finite 
total variation. If the target space X = M, we omit to write it for simplicity. According 
to the Riesz representation Theorem, it is identified to the topological dual of Co{U;X) 
(the space of all continuous functions (p : U ^ X such that {(^ > e} is compact for every 
e > 0), and a weak* topology is dehned according to this duality. The Lebesgue measure 
in R”' is denoted by T”, and the fc-dimensional Hausdorff measure by T-L^. Sometimes, the 
notation ^ will be used instead of for the counting measure, and | • | instead of the 
Lebesgue measure In dimension n, equality of inclusion of sets up to a "H^^^-negligible 
set will be respectively denoted by = and c- 

Given a function u £ L^(C/;R"*) with m > 1. We say that u has an approximate limit 
at X £ [/ if there exists u{x) £ R™ such that 

~ “(®)l = 0- 

JBgix) 

The set Su where this property fails is called the approximate discontinuity set. 
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We say that u has one-sided Lebesgue limits u^{x) e at a; G [/ with respect to a 
direction Vu{x) G := G M”' : |C| = 1} if 

f ^ dy = 0 , 

where B^{x,i'u{x)) := {y G Bg{x) : ±r'u(x) ■ {y — x) > 0}. We will denote by [n](x) := 
u'^{x) — u~{x) the jump of u at x. The jump set Ju of u is defined as the set of points 
X G U such that the one-sided Lebesgue limits with respect to a direction Vuix) exist, and 
in addition tt+(x) / u~{x). Clearly we have Ju C Su- 

2.1. Functions of bounded variation. The space BV{U-,M.'^) of functions of bounded 

variation in U with values in is made of all functions u G such that the dis¬ 
tributional derivative satisfies Du G Wf(C/; The measure Du can be decomposed 

as 

Du = VuC^ + (u+ - u") ® UuB^-^ \-Ju + D%, 
where Vu is the Radon-Nikodym derivative of Du with respect to the Lebesgue measure 
which coincides with the approximate gradient of u. For any 1 < i < m and 1 < j < n, 
we denote by djUi := {Vu)ij the entries of Vu. The measure D'^u is the Cantor part of 
Du which has the property of vanishing on any cr-finite set with respect to the (n — 1)- 
dimensional Hausdorff measure The jump set Ju is a countably ^"■“^-rectifiable 

Borel set, Vu is an approximate unit normal to Ju, and u^{x) are the one-sided Lebesgue 
limits of tt at X G [/ in the direction Vu{x). In addition, we have T-D~^{Su \ Ju) = 0. 

We say that u is a special function of bounded variation, and we write u G SBV (t/; M"*), 
if D^u = 0. If further Vu G for some p > 1, and T-D~^{Ju) < oo, we write 

u G SBVP{U',MD). We refer to [2] for general properties of RR-functions. 

2.2. Functions of bounded deformation. The space BD[U) of functions of bounded 
deformation is made of all vector fields u G L^(f7;]R") whose distributional symmetric 
gradient satisfies 

Dil Du^ 

Eu = -1- e 

This measure can be decomposed as 

(2.1) Eu = e[u)D^ + (n+ - n“) 0 L + E^u. 

In the previons expression, e(u) denotes the absolutely continuous part of Eu with respect 
to For any 1 < < n, we denote by eij{u) = {e{u))ij the entries of e{u). The 

measure E^u is the Cantor part of Eu which has the property to vanish on any cr-hnite 
set with respect to The jump set Ju of u is a countably 'H"“^-rectihable Borel set, 

Vu is an approximate unit normal to Ju, and u^{x) are the one-sided Lebesgue limits of u 
at X G U in the direction Vu{x). If E^u = 0, we say that ■u is a special function of bounded 
deformation and we write u G SBD{U) . We refer to [iQl |35l |38l EHl O [ 8 l [U UHl |T71120] for 
general properties of RH-functions. 

2.3. General conventions. In the sequel we will always work in dimensions 1, 2 or 3. 
Latin indices i, j, k, I, ... (except / and b) take their values in the set {1,2,3} unless 
otherwise indicated. Greek indices a, /3, 7 , ... (except e) take their values in the set {1,2}. 
The repeated index summation convention is systematically used. 
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3. Description of the problem 


3.1. In the original configuration. Let w be a bounded and connected open subset of 
with Lipschitz boundary which denotes the basis of a thin domain occupying the open 
set := u} X (—2e,e) in its reference configuration. We assume that this domain is made 
of the union of a film := cj x (0, e), a bonding layer := cj x [—e, 0], and a substrate 
:= a; X {—2e,—e). Let us underline that the set is not open. Any kinematically 
admissible displacement a : —)• is required to satisfy the boundary condition a = 0 

in In the sequel we shall denote by x' := {xi,X 2 ) the in-plane variable. 

The background behavior of this medium in that of an isotropic linearly elastic material 
whose Lame coefficients are given by 

in Dg. 


The elastic energy associated to a displacement v G satisfying v 

in D= is given by 


(3.1) 



X^eii{v)ejj{v) + 2ifeij{v)eij{v) 


dx. 


0 £^-a.e. 


If the body undergoes cracks, according to the variational approach to fracture (see 
[251 [I3]), the presence of cracks is penalized by means of a surface energy of Griffith type 
where the toughness is given by 


K 


£ 


Kf in Dp 
eub in Df. 


In this case, Sobolev spaces cannot describe admissible displacements since they may 
jump across the cracks. The natural framework is to consider displacements which are 
special functions of bounded deformation. Identifying the cracks with the jump set of the 
displacement, denoted by J„, the surface energy is given by 

(3.2) [ K^dn^. 

The total energy is then given by the sum of the bulk energy, given by m, where e{v) 
is intended as the absolutely continuous part of the strain with respect to the Lebesgue 
measure, and the surface energy, given by (13.2j) . It is well defined for any displacements 
V G SBD{Q,^) satisfying the boundary condition u = 0 C^-a.e. in the substrate D|. 


3.2. In the rescaled configuration. As usual in dimension reduction, we rescale the 
problem on a fixed domain of unit thickness (see m)- We denote by D := D^, Xlf := Dp 
D;, := D^, and D^ := Dj. For every original displacement v G Lr^(D^;M^) (resp. v G 
SBD{Q,‘^)) such that v = 0 C^-a.e. in D=, we define the rescaled displacement u in the 
rescaled configuration by 


Ua{x',X3) = Va{x',eX3), 
U3{x',X3) = eV3{x',SX3) 


for all X = {x',X3) G D. 
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Replacing v by this expression in the energy (j3.ip . and dividing the resulting expression 
by £ yields the following rescaled elastic energy (see m) 


where 

(3.3) Jeiu,nf) := 


dx 


\feaa{u)ei3j3{u) + 2^fea/3{u)eal3{u 

J 2Xfeo,a{u)e33{u) + 4fifea3iu)eo,3iu) 
L (A/ + 2nf)e33{u)e33{u) dx, 


dx 


and 

(3.4) 


J^iu 




+ - 


+ 


Xbeaa{u)epi3{u) + 2fibeal3{u)eap{u 
2Xbeaa{u)e33{u) + Anbea3{u)ea3{u) 

2^2 I (Ab + 2/Xb)e33(u)e33(w) dx. 


Iflb 


dx 

dx 


In the case of cracks, the total energy is obtained by adding the surface energy. In the 
rescaled configuration, it is given by (see [HI EH El E]) 


where 


and 


Ee{u) = Ee{u,nf) -\-Ee{u,ni,), 

“ 1 “ 


j 

1 1 

J 

CO 

/,7unD/ 

\ ^ / 


£/£;('U, ^b) — ^b) E l^b I I ; (^14)3) I 

Ju^^b 


dn^, 


dU\ 


4. Debonding of thin films 

In this section, we assume that the body is purely elastic, i.e., no cracks are allowed. 
Through an asymptotic analysis as the thickness e tends to zero, we rigorously recover a 
reduced two-dimensional model of a thin him system as an elastic membrane on an in¬ 
plane elastic foundation. A similar model has been derived in [341 Theorem 2.1] by means 
of a different method. The original three-dimensional energy Jg : [0,-|-oo] is 

dehned by 

j/\ _ f Je{u,Qf ) + J^{u,Qb) if w G M^) and u = 0 £^-a.e. in 

^ ' I +00 otherwise. 
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while the reduced two dimensional energy Jg : —)• [0, +oo] is given by 


Jo(ii) := < 


_Xf + 2/i/ 


eaa{u)ei3l3{u) + flfeal3{u)eal3{u) 


dx' 


/ IwPdx' 


+00 


if 

otherwise. 


u = (u, 0 ), 
u e Fi(a;;]R 2 )^ 


Our first main result in the following F-convergence type result. 

Theorem 4.1. Let u G then 

• for any sequence (rie)e>o C with ^ u strongly in then 

Jo{u) < liminf Jeiu^); 

£^■0 

• there exists a recovery sequence {u*)e>o C such that u*^ ^ u strongly in 

and 

Jo{u) > limsup j£(u*). 

£—^0 


Proof. Although some parts of the proof are already well known (see [181 Theorem 1.11.2]), 
it will be convenient for us to reproduce the entire argument. 

Step 1. Compactness. Let (ug) C be such that ^ u strongly in 

If liminfg = +oo, there is nothing to prove. We therefore assume that 

liminfg Jeiu^) < oo. Up to a subsequence, there is no loss of generality to suppose that 

^eid^s) — T Jeid^et^b) ^ 


for some constant C > 0 independent of e. The expression (13.3p of the energy in the film 
Llf combined with Korn’s inequality implies that {u^) is actually bounded in H^{Llf,] 


and that ^ u weakly in with u G Contrary to the case of a 

standard linearly elastic plate model (see m), we will show that, thanks to the Dirichlet 
condition in the substrate, the limit displacement u is planar instead of just Kirchhoff- 
Love type. Indeed, using also the expression of the energy ()3.3p - ()3.4p . the fact that = 0 
£3.a.e. in fis, and Poincare’s inequality, we get that 


|(ue) 3 pdx < 




|e33(ue)|^dx < 0, 


so that tt 3 = 0. Thanks again to the bound of the energy in the him (|3.3p . we have 

||ea 3 (w£)||i 2 (f^^) < Ce ^ 0 , 

which shows that ea 3 {u) = 0. It thus follows that d^Ua = —daU^ = 0 which implies that 
Ua{x',X 3 ) = Ua{x') for T^a.e. x G 11/, for some u G We have thus identihed 

the right limit space. 

Step 2. Lower bound. We next derive the lower bound. Up to a further subsequence, 
we may assume that 

£"^633 (Ue) ^ Ca 


--1 


€^0:3 ^ Cq 


weakly in 
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for some functions Ci, ^2 and Cs G Lp‘{^f). Then, by lower semicontinuity of the norm 
with respect to weak convergence, we get that 


lim inf 

>-l 


n f 


[A/(e«a(w) + Cs)^ + ‘2fifeap{u)ea^{u) + 4/r/CaCa + 2/i/C3C3] dx. 


Minimizing with respect to (Ci)C2)C3)) we hnd that the minimal value is attained when 
Ca = 0 and C 3 = -xf^(^aa{u), and thus 


liminf j£(u£, Oj) > 




{u)epp{u) + fifeaj3{u)eafs{ 


u] 


dx. 


A/ + 2/1/ 

We now examine the contribution of the bonding layer. To this aim, according to (I3.4p . 
isolating the only term of order 1 leads to 


Je(li£,n6) > 2/lfe / ea3(li)e«3(li)dx 

J Ui, 

-yJ j j [^3(^^£)2+ '92(u£;)3]dX3 


dx', 


thanks to the Cauchy-Schwarz inequality with respect to the X 3 variable. Since Ug = 0 
£^-a.e. in then 

r-o 

93Ue(x', X3) dx3 = ii£(x', 0 ) for £^-a.e. x' G w, 




where u^{-,0) denotes the trace of on {xs = 0}. On the other hand, setting Ug = 
/°i(«£) 3 (-,a: 3 )dx 3 G H^{u;), we have 

fO 

da{ue) 3 {x',X 3 ) dxs = daU%{x') for £^-a.e. x' G uj. 




Gathering everything, we infer that 


(4.1) Je{Ue,0,b) > 




|(ue)i(x', 0 ) + 5iu|(x')pdx' + ^ / |(u£) 2 (x', 0 ) + ( 92 M|(a:')Pdx' 


According to the trace theorem, and since Ua is independent of X3, we have {us)a{-, 0) —>• Ua 
strongly in Lp'iu). On the other hand, the energy in the bonding layer p.4p together with 
the Cauchy-Schwarz and Poincare inequalities yield 


f |u|pdx'< f |e 33 (iie)p dx < 0 , 

tJ UJ t) 


while p 4 .ll) shows that the sequence (Vug) in bounded in K^). Consequently, Vug ^ 
0 weakly in and combining all the convergences established so far, we deduce 

that 

lim inf (lie, rib) > — f Iwpdx', 

2 

which completes the proof of the lower bound. 
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Step 3. Upper bound. We assume without loss of generality that u = (it, 0) for some 
u G otherwise the limit energy is infinite. We now define a recovery sequence 

For all e > 0, let 

{ {u{x')^£^x^he{x')) ifxGQf, 

(x 3 + l)(it(x'), 0) if a: G 0^, 

0 if X G 

where (/te)£>o is a sequence in C^(lo) such that 


(4.2) 


- ■ X eaaiu) in liin e||V/i£||i 2 (^) = 0. 

Xf + 2flf e^O ^ ’ 


Clearly, tt* G and rt* = 0 C^-a.e. in Qs- Using ()3.3p we have that 

Je{ul,nf) = 1: [ [Xfeaa{u)ei3i3{u) + 2/i/e„/3(it)ea/3(it)] dx 
2Jnf 

+ ^/ [2Xfeaaiu)e‘^hs + fife'^xllVhsl"^] dx 

^ ^ Jq (^/+ 

and according to the convergence properties (14.2p . we get that 

liinj£(ii*,ff/) = ]: [Xfeaa{u)ep^{u) + 2^/e„^(it)e„^(it)] dx' 

e:^0 Z J^ 




Jui + 2 / 1 / 

2 A// 1 / 


2 A? _ _ 1 

^aa (tt)e^^(tt) dx + 


A 2 


2 L A/ + 2 /i/ 


1 


2 Joj . A/ + 2/1/ 

On the other hand, (|3.4I) yields 
-2 


eaa{u)e^l3{u) + 2fifea^{u)ea/3{ 


U 


eaa{u)ej 3 j 3 {u) dx' 


dx'. 


Je(u*, Oft) = ^ / (®3 + 1)^ [AftCaa(n)6/3/3(tt) + 2/tftea/3(■u)eQ ,/3 (u)] dx + ^ / ttatta dx', 

^ ^ Jid 


and thus 


JimJ£«,0/) = ^ J |updx'. 


which completes the proof of the upper bound. 


□ 


5. Transverse cracks in thin films 

In this section, we assume that the body can fracture. We first only address the analysis of 
the film 0/ in order to highlight the appearance of transverse cracks in the reduced model. 
This property is already known in the framework of nonlinear elasticity where energies 
depend on the deformation gradient ig H EH [n]. The difficulty here is to consider a 
linearly elastic material outside the crack so that the energy depends on the elastic strain. 
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5.1. Compactness. From a mathematical point of view, the natural functional setting is 
to consider displacement fields u G For technical reasons, we also assume that 

all the deformations take place in a fixed container K which is a compact subset of 
Therefore, we assume that any displacement is uniformly bounded by some fixed positive 
constant M > 0. 

Throughout this section, we assume that (u£)£>o C SBD{^lf) is a sequence of displace¬ 
ments in the film such that < M, and 

sup < oo. 

£>0 

We establish that any admissible sequence of displacements with uniformly bounded 
energy converges to some limit displacement having a Kirchhoff-Love type structure. 


Proposition 5.1. Up to a subsequence, there exists u G SBD{Qf) such 

that 

i) Us ^ u strongly in and ^ u weakly* in 

a) e{us) e{u) weakly in 

Hi) ea 3 {u) = 633 ( 1 *) = 0 C^-a.e. in Qf and {i'u )3 = 0 T-L^-a.e. on 


Proof. From the hypotheses and the definition of Es{-,^f)j have that 
ll’*^e||l/°°( 0 /) T ll®(’^e) + PL {Jue Fl 11/) ^ C, 

for some constant C > 0 independent of e. According to the compactness theorem in 
SBD [ini Theorem 1.1], we deduce the existence of a subsequence (not relabeled) and a 
function u G SBD{^}f) such that ^ u strongly in L^(n/;M^), ^ u weakly* in 

e(us) e{u) weakly in L^(ll/; Mgy^)]), and 


n'^iJuunf) 


(5.1) 


€—^0 

< liminf 
£—^0 




Using the expression of the energy in the him, we deduce that 


dn‘^. 


(5-2) \\ea3{ue)\\L\nf)+ [ |(z"«e) 3 |d'H^ < Ce 

and 

l|e33('“e)llL2(Qy) < Ce^ 

for some C independent of e. Using the lower semicontinuity of the left hand side of both 
equations with respect to the convergences established for (i*£)e>o (see [TOl Corollary 1 . 2 ]) 
we conclude that 603 ( 1 *) = esslu) = 0 £^-a.e. in 11/, and that (nu )3 = 0 PL^-a.e. on 
JuUiQf. □ 


In the sequel, u denotes a displacement as in the conclusion of Proposition 15.11 Our 
next goal is to get a more precise structure of such displacements. Contrary to the case 
of linear elasticity (see [18]) or linearly elastic-perfectly plastic plates (see [ 22 ]), they in 
general are not of Kirchhoff-Love type (i.e. such that Ei^u = 0) since we do not control 
the full distributional strain Eu. In particular, the singular part of the shearing strain 
Ea 3 U is given by which might not vanish. However, we shall prove below 
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that they have the same structure in the sense that the transverse displacement only 
depends on the planar variable x', while the in-plane displacement {ui,U 2 ) is affine with 
respect to the transverse variable x^. 


Proposition 5.2. Let u £ SBD{Q,f) n be such that ei 3 (u) = 0 C^-a.e. in 

Llf, and {i'u )3 = 0 T-L^-a.e. on Then the following properties hold: 

• the function is independent ofx^ and it (is identified to a function which) belongs 
to SBV{uj) n L°°{u}). In addition, its approximate gradient Vus = (5i 1 x 3 , 521 x 3 ) £ 

• for C^-a.e. {x',xs) £ 


(5.3) 


Ua{x',X 3 ) 


Ua{x') + 




daU^ix'), 


where Ua ■= Jq Ua(-, X 3 ) dxs, and u := {ui,U 2 ) E SBD{u}) H 
• Ju = {Ju u Ju3 u Jvug) X ( 0 , 1 ); 


Proof. Step 1. First of all, by virtue of (12.11) . the distributional derivative of U 3 with 
respect to X 3 satisfies 

D3U3 = E33U = e33{u)C^ + [lx]3(r'u)3^^L Ju = 0 . 

This implies that 1 x 3 is independent of X3, and that it can be identified to a function defined 
on w. 


Step 2. We next show that 1 x 3 £ SBV{co) and that formula (15.31) holds. This will be 
obtained thanks to a suitable mollification of u. We hrst extend u to the whole space in 
the following way: since the trace of an SBD{Q,f) function belongs to (see [H 

Theorem 3.2]), according to Gagliardo’s Theorem, u may be extended to by a function, 

still denoted by u, that is compactly supported in and such that u £ M^) 

with|£u|(an;) = o. ^ \ i 

Let X £ be an even and non negative function such that J]gx(^)df = 1 and 

Suppx C (—1,1). For all x = {xi,X 2 ,X 3 ) = (x',X 3 ) £ M^, we define p(x') := x(xi)x(x 2 ) 
and p(x) := x(xi)x(x 2 )x(x 3 )- We then denote by psix') = 5~'^p{x'/5) a sequence of 
two-dimensional mollihers, and by p 5 {x) = 6 ~^p{x/ 6 ) a sequence of three-dimensional 
mollihers. Since u* ps ^ C^(M^;M^) and 

53(11 * ps)a = 2 ea 3 (lX * ps) - da{u * ps)3, 
it follows from the fundamental theorem of calculus that for each [x',X 3 ) £ Llj, 


(5.4) Ua* psix' ,X3) = Ua* psix',0)+ 2 e „3 (ix * p^) (x', s) ds 

Jo 

- / da{u3 * P5)(x',s)ds. 
Jo 

Let us study each of the above terms separately. The term in the left hand side of (|5.4p 
clearly satisfies Ua* ps ^ Ua strongly in L^(flp), and thus (for a suitable subsequence) 


(5.5) Ua*ps^Ua C^-a.e. in Llf. 

Concerning the first term on the right-hand side of ()5.4p . standard properties of convo¬ 
lution of measures ensure that Eus ^ Eu weakly* in A4(M^) and ji^ii^KM^) —)• |£'ix|(M^). 
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Therefore, since \Eu\{dQf) = 0, we deduce that \Eu^\{nf) —)• \Eu\{il.f) which implies, by 
the continuity property of the trace (see [SI Proposition 3.4]) that Ua* ps ^ Ua strongly 
in Thus, denoting by ii+(-,0) the upper trace of Uq, on w x {0}, there is a 

subsequence such that 

(5.6) Ua * ps{',0) ^ E^-a..e. in w. 

Regarding the second term on the right-hand side of (15.41) . we have Casiu * ps) = 
(EaSU) * PS with EaSU = 

£{x',X 3 ) := / Casiu* Ps){x',s)ds 

Jo 

1 r 

= a / Psix'-y',s-y3)[u]3{y){uu)a{y)d'H^{y)ds. 

Jo J Ju 

Since u G with < M, then |[ii]| < 2M which leads to 


|^^(a;',a;3)| < M 

= M 



r-l 

'0 JJr, 

. 

'.J-u. Jo 


Psix' -y',s- ys) dn‘^{y) ds 
Psix' -y',s- y3)dsdn^iy), 


where we used Fubini’s Theorem in the last equality. We next denote by Q'ix',5) := 
x' + (—6, 5)^ the open square of (parallel to the coordinate axis) centered at x' and 
of edge length 2S. Observing that psix' — y',s — ys) = 0 if y' 0 Q'ix',6) and that 
Ps{x' -y',s- ys) = Psix' - y')6~^xiis - yo)/^) with xiO dt = 1, we get that 

\£ix',X 3 )\ <M f psix' -y') ( [ (5-^x((s-2/3)/<5)ds) d^^(y) 

J J-u.n[Q'ix',5)x{o,i)] \Jm J 


= M 


>1 

J.h. 


Psix' - y')d'}i^iy). 


' .J-u.n[Q'{x',S)x{0,l)] 

For any Borel set B C co, let us define the measure piB) := JJ?iJu 0 (R x (0,1))) which is 
nothing but the push-forward of Ju by the orthogonal projection vr : —)• x { 0 }. 

Note that p is concentrated on 7 r(Ju) since //(a;\ 7 r(J^)) = 'H‘^iJu^[ioj\'^iJu)) x (0,1)]) = 0. 
On the other hand, the generalized coarea formula (see [21 Theorem 293]) yields 


C\7:iJu)) < 


^“(J^nvr ^ix'))dx'= / \ivu)‘i\dJE = d. 


Ju 


’ 'e{Ju) 

Therefore, p and are mutually singular which ensures that the Radon-Nikodym deriv- 

dji 

dZJ 

|2 


ative 4 ^ix') = 0 at £^-a.e. x' G oj. R follows that for T^-a.e. x' G ca, 


sup |T(a:',a: 3 )| < ^ 

X3e(0,i) 


(52 


and thus, in particular. 


(5.7) 


rX3 

Jo ' 


e-aoiu * Ps)ix'■, s) ds —)• 0 for £^-a.e. ix', 3:3) G fi/. 
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For what concerns the last term on the right-hand side of (I5.4p . since is independent 
of X 3 , we infer that U 3 * ps is independent of X 3 as well since U 3 * ps{x) = U 3 * psix') for 
all X £ Therefore, 

da{u3 * P5){x', s) ds = X 3 da{u 3 * Ps){x'), 
and (j5.4l) - (15.7p thus imply that 



da{u3 * P5){x^) 


^I+(x^0) - Ua{x\x 3 ) 

X 3 


ijjaix') for C^-a.e. {x',X 3 ) £ ilf. 


That V'a only depends on x' is due to the fact that the left-hand side only depends on x'. 
Moreover, since u;^(-,0) £ and Uq,(-,X 3 ) £ L?‘{oj) for a.e. X 3 £ (0,1), we deduce that 

From the last formula we get that 


(5.8) Ua(x',X 3 ) =tt+(x',0) - X 3 'ljjaix'), 

which in particular implies that D 3 Ua = and 

DaU3 — D3'Ua T — i/’q./I “1“ \'^3{yu)aH I— Jw 


As a consequence, the distributional derivative in Vlf of 1 x 3 is a bounded Radon measure 
in flj, and therefore U 3 £ BV{Qf). Since the singular part of the above measure is con¬ 
centrated on Ju which is fj-hnite with respect to we deduce thanks to [21 Proposition 
3.92] that U 3 £ Finally, since U 3 is independent of X 3 , we actually infer that 

U 3 £ SBV{uj). In addition, by uniqueness of the Lebesgue decomposition, it follows that 

Tpci — daU 3 , ['fx\ 3 (^l'u')a'H L Ju — [ns](^113)o'^ ^ [“Ajs ^ 1)] 


SO that 


(5.9) Ju 3 X {0,'^) C Ju- 

Integrating relation (15.8p with respect to X 3 yields 


u. 


I 

t{x') := / Ua{x', X3 ) dx3 = Uq(x', 0 ) — -doiU 3 {x') for £^-a.e. x' £ cu, 
Jo 2 


from where (15.31) follows. 

Step 3. Let us prove that the approximate gradient of U 3 , denoted by VU 3 := 
(9i 1 x 3 , 921 x 3 ), and the averaged planar displacement u := (xxi,xx 2 ) belong to BD{u}). For 


|r2x2 
^sym 

orem 3.2]), we infer that 


any p £ C^(a;; Mgyj^), according to the integration by parts formula in BD (see | 8 l The- 


- / dpifapUadx' = - / dpPapUadx = / / <y9„^XX„U^ d'H^. 

J u) ^ f / 

Since (^ = 0 in a neighborhood of dio x (0,1) and u = ±63 on to x {0,1}, we get that the 
boundary term in the previous expression is zero. Therefore 

(5.10) - / dpPapUocdx' = / Papeap{u)dx+ / Pa.pi\u\Ql'u)apd'H^ 

«/ tjJ J J u 

which shows that u £ BD{ui). According to slicing properties of BD functions (see [U 
Proposition 3.4]), for £^-a.e. X 3 £ (0,1), the function (xxi(-,X 3 ),XX 2 (-,X 3 )) £ BD{uj) so 
that relation (15.3p yields in turn that Vxx 3 £ BD{io). 
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Step 4. We next establish that J-^ = (JitU J^gU x (0,1). To this aim, let us dehne 
the functions v := {ui,U 2 ,u^) and g := ( 9 i 1x3, 921*3) 0). Since 113 G SBV{uj), u G BD{uj) 
and Vu 3 G BDioj), then clearly both v, g ^ BD{Qf), and 


(5.11) Jg = Jv.g X (0,1). 

Moreover [21 Proposition 3.92 (b)] and [H Proposition 3.5] imply that 


T- — 

^11. - 


x' G io : lim sup 


\Eu\m^')) 


> 0 


= < x' G oj : lim sup 
£>^■0 


\Dus\mx')) 


> 0 


and 


^ 0 „ X G 11 




lim sup 


\Ev\iB,ix)) 


> 0 


where B'^{x') stands for the two-dimensional open ball of center x' and radius g, while 
Bg{x) stands for the three-dimensional open ball of center x and radius g. Since v is 
independent of X3, then 


(5.12) J^^(J^UJ„g)x(0,l). 

According to (j5.11|l and (j5.12|l . it is thus enough to show that = J^, U Jg. 
Let us also dehne the sets 

I Q^O 

©g := |x G : lim sup \^ 3 \{Be{x)) ^ ^ 

I (?^0 Q 


and recall that, according again to [H Proposition 3.5], ©„ = and Qg = Jg. Using the 
expression of the displacement ()5.3p . we have u = v + {\ — x^)g. Since u G L°°(11;M^), 
then V G L°°(a;;M^) as well, and the previous relation yields g G L°°(a;;M^) with 

If 

lim \g\dy = 0 for all x G Hf. 

r JBgix) 

Consequently since Eu = Ev + {^—X‘i)Eg — e^Qg, we deduce that ll/\(0„U0g) C Qf\Qu 
i.e. &u C 0.U U 0g and 

(5.13) JuCJv'J Jg- 

We now prove the converse inclusion. From the relations v = u+{x 3 — ^)g and (^ —X 3 )fir = 
u — V, and the fact that g is independent of X3, we similarly obtain that 0.y C ©^ U 0g 
and 0g C Qu U Qv which imply that 

(5.14) Jv\JgCJu, Jg\JvCJu. 


It thus remains to prove that 


(5.15) Jv(^JgC Ju- 
According to (15.91) . (|5.11l) and (I5.12p . we get 

(J„ n Jg) \ ^ {[{Ju n Jv„g) X (0,1)] \Juc [{Ju n Jy^) \ Su^] x (o, i), 

where we used that, since U 3 G SBV{uj), then J^g = Su^. Assume by contradiction that 

(5.16) ^^{{Jv n Jg) \ J„) > 0, 
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then there is some x = {x',x^) G (J^ H Jg) \ Ju with x' G (J^ H \ Su^ such that 

^uix') = ±r'vu3(x'). Let us assume without loss of generality that i'u(x') = vs/usix') =: 
i^{x'), the other case can be dealt with similarly. Since x' is a Lebesgue point of 1x3, then 
the one-sided Lebesgue limits of U3 at x' in the direction u(x') are equal and coincide 
with its approximate limit. On the other hand, since x' G Ja H Jvua; then the functions 
u and Vu3 admit one-sided Lebesgue limits at x' in the direction i'(x'). Next, from the 
expression ()5.3I) of the displacement, we deduce that for all a G {1,2}, the functions Uq, 
admit as well one-sided Lebesgue limits at x in the direction (u(x'),0). Gathering all 
previous informations, we get that the full displacement u admits one-sided Lebesgue 
limits at x in the direction (u(x'),0). Using the fact that x 0 Ju, we infer that necessarily 
[it](x) = 0, and thus, using again ()5.3p yields 

(5.17) + Q - [Sails](ic') = 0 for all a G {1,2}. 

We observe that, by (15.lip and (I5.12p . the sets and Jg are invariant in the transverse 

direction, and consequently {x', 1/3) G Jy O Jg for any G (0,1). Therefore if (x', 2/3) 0 Ju 

for some 1/3 7^ X3, then reproducing the same argument than above implies that 

[Ua]{x') + Q “ [Sails](x') = 0 for all a G {1,2}. 

Subtracting the previous relation to (I5.17P yields [ua]{x') = [daU^\{x') = 0 for all a G 
{1,2}, which is against the fact that x' G J^ n Jvua- As a consequence, (x',?/3) G Ju for 
all 1/3 G (0,1) with 1/3 7^ X3. In addition, since x' G Jy^a, there is some a G {1,2} such 
that [Soils] (x') 7^ 0, and X3 is therefore given by 

1 [Ua\{x') 

^ 2 [Solis] (x')’ 

Consequently, we have proved that 

{Jv^Jg)\Ju C y I {x',x^) : x' G Jii n Jy„3, [Solis](a;') 7^ 0, X3 = ^ ) 1 =. 

I ^ [DalisH.^; j J 

The set A is Borel measurable, and, for each x' G Ju n Jy^g, its transverse section passing 
through x', denoted by := {xs G (0,1) : (x'jXs) G A} is reduced to at most two 
points. Since the set Ju H Jy^a is countably 77^-rectifiable, [231 Theorem 3.2.23] ensures 
that Ti.^ L ((Ju n Jyua) X (0,1)) = ("H^ L (Ja n Jsjui)) ® (71^ L (0,1)), and Fubini’s Theorem 
yields 

772(A) = f C\A^')dn\x') = 0 , 

which is against (I5.16p . and therefore completes the proof of (I5.15p . Gathering (I5.13P - 
()5.15l) leads to Ju = J^ U Jg, and thus Ju = (Ja U J^a U Jsjuz) ^ (0, !)• 


Step 5. We complete the proof of the proposition by establishing that u and Vits are 
actually SBD{ui) functions. Indeed, since we know that J^^ = T x (0,1) for some countably 
77^-rectifiable set T C w, equation (15.1011 reads 

- / d^ipa^Uadx' = / (fapi Ca^iu) dx^] dx'-7 (fa/si (M © i"r)a/i dx3 ) d77^ 

J oj J u) \Jo / Jt \Jo j 
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which implies that Ca/siu) = eQ,^(u)(-, X 3 ) dxs by uniqueness of the Lebesgue decompo¬ 

sition, and that the singular part of Eu is concentrated on a countably ^^-rectifiable set. 
It follows from [TJ Proposition 4.7] that u G SBD{u) and the same can be said, therefore, 
first for Vtts and then for («]*'(•, 0 ), u^(-, 0 )). □ 

Propositions 15.11 and 15.21 suggest one to define the limiting space of all kinematically 
admissible displacements by 

Akl ■■= l^u£SBD{nf):\\u\\Loo(^^^-^<M, U3£SBV{uj)nL°°{uj) 
with Vu 3 e SBD{uj) n L°°(a;;M2), 

(5.18) Ua{x', X3) = Ua(x') + (2 ~ X‘^daU 3 {x') for £^-a.e. x = {x', X 3 ) G fl/, 

where u := {ui,U 2 ) G SBD{u) H M^), 

(5.19) and Ju = (Jn U J „3 U Jvug) X (0,1) >. 


5.2. P-limit in the film. For each e > 0, let us define the functionals and : 
L^(nj;M^) [ 0 ,-|-oo] by 


and 


if 11 G S’.BZl(fl/) and ||u||£,oo(Q^) 
^ ^ ^ ' I +00 otherwise. 


< M, 


^o(^) := ^ 


_Xf +2nf 


eaa{u)epi3{u) + fj.fea/3{u)eal3{u) 


+1^ 




dx' 


_Xf + 2/ij 
d-K/B {Ju U Ju 3 U 
-boo 


eaa(V'U3)e/3/3(VM3) -b /i/ea/3(Vu3)e„,g(Vu3) 


dx' 


if It G Akl, 
otherwise. 


Theorem 5.1. The sequence of functionals (£^/)£>o T-converges to £q with respect to the 
strong L'^{Qf,'R^)-topology. 

Proof Step 1. We start by deriving a lower bound inequality, i.e., for any u G L‘^{Qf] M^) 
and any sequence C L^(fIy;M^) such that ^ u strongly in L^(nj;M^), then 

lim inf {u^) > £q (it). 

r 

If liminfefie (it^) = -boo, the result is obvious. Otherwise, up to a subsequence, we can 
assume that 

lim £^ (itj) = lim inf £^ (it^) < 00 . 

£—>■0 e^O 

By virtue of the above energy bound, we can assume without loss of generality that the 
conclusions of Propositions 15.11 and 15.21 hold so that it G Akl- Using a very similar 









REDUCED MODELS FOR LINEARLY ELASTIC THIN FILMS 


19 


argument than that used in the proof of the lower bound in Theorem 14.11 combined with 
the lower semicontinuity of the surface energy established in ()5.1j) . we obtain that 

According to (|5.18p . we get that 


liminf flj) > 


_Xf + 2nf 


eaa{u)ei3p{u) + Hfeal3{u)eal3{u) 


/ eai3iu)eai3{u)dx = eai3{u)eap{u) 

J J Q,f 


+ 2 { - - Xa ) ea/3(w)e„^(Vu3) 


+ (2 “^3] ea/3(Vu3)e„/3(Vu3) 


dx 


= / eai3{u)eai3{u) d®' + — / e„y3(Vu3)e„^(Vu3)dx', 


and similarly for the other term 

f eaa{u)ep0{u) dx = I eaa{u)ei3i3{u) dx'+ ^ f eaa(Vtt3)e/3/3(Vu3)dx'. 

Ju) jUJ 

Therefore ()5.19p yields the announced energy lower bound. 


Step 2. We next derive an upper bound through the construction of a recovery 
sequence, he., for every u G there exists a recovery sequence {u*)^yQ C 

such that u* ^ u strongly in and 

limsupT/(u*) < £q{u). 

£—^0 


If u 0 Akl, then £q{u) = +00 and the result is obvious. It therefore suffices to assume 
that u G Akl- We now define a recovery sequence (w*)£>o. For C^-a.e. x = (x', X 3 ) G ri/ 
and all e > 0 , let 

u*(x',X 3 ) = c^{u{x) + { 0 , 0 ,e‘^X 3 he{x'))), 
where {he)e>o is a sequence in C^(a;) such that 


(5.20) hs -- —^2-— Baaiu) in lini e||V/i£||i 2 (^) = lim e\\he\\L^(^^) = 0, 

e->-0 ^ ^ ' 

and Ce := Mj{M + £^\\he\\L°°{^^))- Clearly, u* G SBD{£lf) and < M. Using 

P3.3I) we get that 


= % / [Xfeaa{u)epji{u)+ 2^ifec,p{u)ea0{u)\dx 
^ JUf 

(? f 

+ ^y [2\feaa{u)£^he + ^lfe^xl\Vhe\^] dx 
Thus, since —)• 1 and according to the convergence properties (|5.20l) . we get that 


lim 



Xf + 2fif 


eaa{u)epii{u) + 2flfeal3{u)eal3{u) 


dx. 
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Concerning the surface energy, since = Ju 
that 


/ 



V e / 


dn^ 


which completes the proof of the upper bound. 


= {Ju U Ju 3 U Jvms) X (0,1) it follows 
= H^iJu U Jus U Jvus), 

□ 


6. Multifissuration: debonding and delamination vs transverse cracks 

In this section, we consider the full model of a film f deposited on a substrate Qg through 
a bonding layer fife, and we assume that both fij and fib can crack. 

6.1. The anti-plane case. Following [32] . it is assumed that the geometry is invariant 
in the direction 62 , i.e., uj = I x M, where / is a bounded open interval, and that the 
admissible displacements take the form 

u{x) = u{xi,X3)e2. 

In this case the elastic energy reduces to 

Je{u) = ^ [ + e"^|53ttp) dxi dx3 

^ i/x(0,I) 

+ ^[ (e^|<9in|^ + |53ttp)dxidx3, 

^ J/x(-l, 0 ) 

and the total energy is given by 

Es{u) := Je{u) + Kf f \{{i^u)i,£~^{i^u) 3)\ dn^ 

./J„n[/x(o,i)] 

+ Kb [ |(e(t'n)l, (z^n)3)| 

./j„n[/x[-i,o]] 

The natural functional setting is to consider (scalar) displacements in the class 

A := {u £ SBV{I X (—2,1)) : u = 0 £^-a.e. in / x (—2, —1) and ||u||j;,oo(/x(o,i)) < M}, 

where M > 0 is an arbitrary fixed constant. 

In [32], the following one-dimensional energy, defined for all u G SBV{I), was proposed 
as an approximation of the previous two-dimensional energy 

Eo{u) ^ f \u'\‘^dxi + ^ f \u\'^ dxi + Kfi^{Ju) + KbE^{Au), 

^ Jl ^ Jl\Au 

where := {|tt| > ^^/2Kb/Jj^} is the delamination set. An easy adaptation of the proof of 
[331 Theorem A.l] justifies rigorously this conjecture through the following F-convergence 
type result. 

Theorem 6.1. Let u G SBV{I), then 

• for any sequence (u£)£>o C A satisfying ^ u strongly in LJ‘{I x (0,1)), then 

Eq{u) < liminf .&£(«£:); 

£^>0 
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• there exists a recovery sequence (n*)e>o C A such that u* ^ u strongly in L‘^{I x 
(0,1)), and 


Eo{u) > liminf (m*). 
£^0 


Let us observe that if is a sequence of minimizers of (under suitable loadings), the 
(characteristic function of the) delamination set is constructed as the L^-limit of the 
orthogonal projection of the jump sets Ju^ onto the mid-surface { 3:3 = 0}. In particular, 
the vertical cracks in the bonding layer do not contribute to delamination. 


6.2. The general case. We conjecture that Theorem 16.11 can be extended to the general 
three-dimensional vectorial case. In this situation, the space of kinematically admissible 
displacements is given by 

A := £ SBD{il.) : u = 0 C^-a.e. on Qg, and < m|. 

Let us define the energy functionals Eg and Eq : LP‘{EL; M^) [0, -|-oo] by 

Eg{u) u ^ A, 


( 6 . 1 ) 

and 


Ee{u) := 


-|-oo otherwise. 


Eq{u) := < 


_Xf -h 2/i/ 


{u)ei3g{u) -F yLfeaii{u)eap{ 


u 


+ 

+ 

-|-oo 


12 


^ Ju}\A 


_Xf + 2fif 


dx' 


e«a(Vir3)e^/3(Vu3) -F ///e„^(Vu3)e„y3(V«3) 


dx' 


|wP dx' -F U J„3 U Jvus) + KbC^{A) if n e Arl, 


otherwise, 


where the delamination set is defined by 
( 6 . 2 ) 


A := |x' £ CO : |■u(x')| > u {x' £ cj : U3 / 0}. 

We expect Eq to be the L-limit of as e —)• 0, but have been unable to prove the 
corresponding lower bound inequality: 


Conjecture 6.1. If u £ L^(fl;M^) and {ug )g>o C L2(L!; 
strongly to u in L^(r2j;M^), then 

Eo(u) < liminf ^^(ug). 
£^■0 


IS any sequence converging 


Our aim here is only to prove the L-limsup inequality and to present some partial results 
and techniques which could be relevant in future investigations of this problem. 

Proposition 6.1. For every u £ L^(n;M^), there exists a sequence {u* )g>0 C L2(L?;M3) 
such that u* ^ u strongly in L^(nj;M^), and 

Eo{u) > lim sup £^£(u*). 

£—^0 
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Proof. If u 0 AkLj then So{u) = +oo and there is nothing to prove. Therefore, we assume 
from now on that u G Akl- 

Step 1. In order to define the recovery sequence, we need several approximation steps. 
We start by approximating the delamination set defined in (16.2p by a sequence of sets 
of finite perimeter. Let {pm)meN be a standard sequence of mollifiers in and set 
Xm ■= Pm*XA- We know that Xm —^ XA strongly in L^(a;). Set 

■— \J\\Xm ~ XAllLl(tj) ^ 0. 

By the coarea formula [21 Theorem 3.40], for every m G N large enough, there exists 
\ <1 — 5m such that 

•— {3^ S U1 . Xmipj ) ^ im} 


has finite perimeter. We claim that 
(6.3) XA^XA in L^(a;). 


Indeed, 


and 


\ A) < ^ f Xm{x') dx' < ^ 
tm J Am\A 



Xa \ dx' —^ 0, 


£2(A \ Am) < C‘^i{x' G A : xa(x') = 1 and Xm(x') < 1 - 5m}) 

< ]- [ Ixa(x') - Xm(x')| dx'< 5m ^ 0, 

Om J A 

hence jjxA - XAm\\p{uj) = \ A) + ^^(A \ Am) —^ 0. In addition, it is possible 

to find a sequence £m g such that emP^id*Am.) With a slight abuse of 

notation, we refer to the sequences (Sm) and (Am) simply as (e) and (A^) and henceforth 
assume that 


(6.4) lim£nUd*AA = 0. 

£^•0 

We next approximate the displacement u. Indeed, according to m Theorem 3] (see also 
[3U1 Theorem 3]), there exists a sequence G SBV (cj; M^) such that —)■ u strongly 

in ^ strongly in 4 .^( 0 ;; H^iJue \ Ju) + B^iJu \ Jue) 0> 

and lluellLoo^^) < ll'ujl^oo^^^). Let us define for a.e. x' G w and all X 3 G (0,1), 

{Ue)a{x',X3) := (Ue)a(x') + (^ j OaUsix'), {Ue)3{x',X3) := U^ix') 


SO that Ue G SBD{Q,f), and ||ue||ioo(Q^) < < M. 

As in the proof of Theorem 15.11 we consider a sequence {hs)e>o C C)T(w) satisfying 

(CTD . 

We now define the recovery sequence by setting, for all e > 0 and for £^-a.e. x = 
(x',X 3 ) G 0, 


U*(x',X3) 


Ce{ue{x) + (0,0,e^X3/i£(x')) if (x',X 3 ) G ri/, 

< C£(x 3 + l)(we(x'),0) if (x',X 3 ) G (a; \ AJ X [-1,0], 

^0 if (x',X3) G (Ag X [-1,0]) U 
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where ) ' Since the set has finite perimeter in u and S SBV{ui; M^), 

then u^Xuj\Ae ^ < 55 ^( 0 ;; M^), and thus u* S SBD{Q,) and ti* = 0 C^-a.e. in In addi¬ 
tion, the fact that \\us\\L<=‘>{nf) ^ ^ yields ||ite||L°o(o^) < M as well so that u* G A. The 
sequence (u*)£>o is thus admissible, and clearly u* ^ u strongly in L‘^{Qf]'K^). 

Step 2. Using the convergence properties of tie, a similar argument than in the proof 
of Theorem O leads to 

limsupiile(w*, flj) <£ q { u ). 

£—)'0 

It thus remains to compute the energy associated to this sequence in the bonding layer. 
First, the bulk energy in the bonding layer gives 


Jeiu*,nb) = 


cy 


f 

J (l 


(a;\Ae)x(-l,0) 


(X3 + 1)^ Xbeaa{Ue)e^i3{Ue) + 2/Xf,e«/3(Ue)e«/3(u£ 


dx 


+ 


CefJ-b 


f |w£p dx' —^ / 

J Lj\Ag 2 J u)\A 


IllP dx'. 


Concerning the surface energy in the bonding layer, we first observe that for each e > 0, 


Tit* n £ib c 


J^,x[-1,0] U A,x{0} U ({(iX3,Vir3) /0}\A,)x{0} 


U 




A,x[-1,0] 


where d*Ag stands for the reduced boundary of A;. [21 Definition 3.54]. Let us observe 
that cj \ A C {u 3 = 0} c VU 3 ) = 0} since, by locality of the approximate gradient, 
Vn 3 = 0 £^-a.e. in {rt 3 = 0} (see jH Proposition 3.73 (c)]). Then 


limsup [ \ 

£^>0 7,/^*n06 


< limsup 
£—^0 


en\Ju,) + £2(A,) + £2({(^3^ V 1 X 3 ) / 0} \ A,) + en\d*A,) 


thanks to (16.21) . (16.31) and (16.41) . 


£2(A), 

□ 


6.2.1. Partial results for the lower bound. Let u G and iug),>o c L'^iny) 

be a sequence such that Ug ^ u strongly in If liminf^ = -|-oo there is 

nothing to prove. Otherwise by (16.ip . up to a subsequence, we can assume without loss of 
generality that £ A, and that 

(6.5) sup£(£(iie) < -boo. 

£>0 

As a consequence, all the compactness results in the film D f established in section 15.11 
hold. In particular, Propositions 15.1] and [A2] show that u G Akl, and the lower bound 
established in Theorem 15.11 yields the terms in £o{u) corresponding to the energy in £lf. 
The main problem is to deal with the bonding layer. Following the scalar case treated 
in |33] . it is enough to show that the energy in fib is bounded from below by some func¬ 
tional where the delamination set is replaced by a function 6 G [0,1]), which can 

be interpreted as a delamination volume fraction density. On {9 = 1}, the him is entirely 
debonded from the substrate, while on {9 = 0} it continuously accommodates the pre¬ 
scribed zero displacement on the substrate exactly as in the Sobolev case fTheorem 14.11) . 
All intermediate states are contained in the set {0 < 0 < 1}. 
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Proposition 6.2. Assume there exists 6 G (w; [0,1]) such that (1 — 9)u3 = 0 C^-a.e. 
in oj, and 


( 6 . 6 ) 

Then 


2 


2 


(1 — dx'+ Kfe / < liminf Ofe). 

■up dx'+ Kb£^(A) < liminf £^£(u£, rife), 


> lj\A 


£—^0 


where A be the delamination set defined in 
Proof. By assumption, we have that 

^(1 - r/)|n(x')|^ + Kbr]) dx' < liminf E^{ue, Ofe). 

The result follows by solving the pointwise minimization problem explicitly. 


/ min 

Juj {»7e[0,l]:(l-r;)u3(x')=0} 


□ 


The main point is to construct such a function 9. As in the scalar case [33], 9 is supposed 
to be obtained as the L°°(ti;)-weak* limit of a sequence (yAe)£>o of suitable measurable 
sets Ag C ui. However, it is unclear what is the right notion of an e-delamination set A^ 
in the vectorial case. In particular, the following example shows that vertical cracks in 
the bonding layer cannot be neglected, so it is not enough to define A;, as the orthogonal 
projection of Ju^ onto the mid-plane uj x {0}, as in the anti-plane and in the Sobolev case 
(Thm. [6Tl [33] Prop. B.2], and Thm l4.ip . 


some 


Example 6.3 (Microstructure example). Suppose that u = (0,1)^ and e = ^ for 
A G N. In the film, set 

Ui;{x) = u{x) = (0,^,0) for all x G Hj. 

In Hfe set, for each i = 0,..., A — 1 and all 2ie < X 2 < (2i -|- 2)e, —1 < X 3 < 0 < xi < 1, 

' X2 — 2is 


Ueixi,X2,X3) = 0 ,i{l + xfij,iev 


1 + 3:3 


where v G ^(0, 2) x (0,1)^ is any function such that v{s, 0) = v{s, 1) = 0 Vs G [0,1] and 

q ■■='f [ {{1 + dsv)'^ + 2dtv‘^) dsdt < 1. 

J s=0 J t=0 


f s=0 J t=0 

If Ag is defined as '^{Jue H H); then 


(6.7) / (2;Ufee«3(u£)ea3(rie) + e Hbe33{ue)e33{us)) dx -T KbC‘^{A^) 

J Qh 


q^ibf 


On the other hand, if A = {|ti| > ^/2Kb/fib} is the expected limit delamination set, then 


( 6 . 8 ) 


/ y 

Juj\A 2. 


UaUa dx' + Kfe£^(A) = 


Kb 


if e< 
if 


Choosing £ G ( \ shows that (j6.8l) would not be a lower bound for (j6.7p . □ 




qi^bj 
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Regardless of the notion of an e-delamination set one tries to define, it is convenient 
to impose that it should contain the set 

Pe ■= n f^/), 

where tt : —)• 7t(x) := x', is the orthogonal projection onto R^ x {0}. On the one 

hand, there is no loss of generality in doing this, since it converges to a Lebesgue negligible 
set. Indeed, according to the coarea formula (see [21 Theorem 2.93]) and the surface energy 
bound (|6.5p in the film, we have 

C\Pe)<l 'H^{Juennfn7r-\x'))dx'= f Ku^JsIchR^ < o. 

On the other hand, excluding Pi; enables one to slightly improve the convergences in the 
film, as in the following lemma which proves the convergence of the planar gradient of the 
anti-plane displacement. 

It will be assumed henceforth that U; G and that Ju^ is closed in Q and 

contained in a hnite union of closed connected pieces of hyper surfaces. In doing this 
no generality is lost, thanks to the density result in SBD of m Thm. I]. In particular, 
we have that U; G \ P;) x (0,1);R^). 


Lemma 6.4. Let (A£)£>o he a sequence of closed sets he such that P; C Ag for each 
e > 0. Assume that there exists a function 6 G [0,1]) such that XAe ^ ^ weakly* in 

L°°{uj), and (1 — 9 )u 3 = 0 C^-a.e. in ui. Then 

Xuj\A,da{ue )3 ^ 0 Weakly* in L^{u; 1)). 

Proof. First note that for C^-a.e. x' ^ P; and £^-a.e. X 3 G (0,1) 

r^3 

:= / (9„('W£)3(a;',s)ds-F (u£)„(x) - (u£)+(x',0) 

Jo 


= 2 


ea3(ii£)(3:',s)ds. 


= / [daiUe)3{x',s)+d3{Ue)o.ix',s)]ds 

Jo 

Thanks to the bulk energy bound (|6.5p in the him (see also (|5.2p i. we have that 
(6-9) llCQllL 2 ((^\Ae)x(o,i)) < 2 ||e„ 3 (u£)||i 2 (Q^) < Ce ^ 0. 

Integrating (16.91) we obtain that also ||CallL 2 (tj\Ae) 0; where 

^1 /-Xs 

Caix') ■= CUx', X3) dX3 = / da{Ue)3ix',s)dsdX3 + {U;)aix') - {U;)'^{x',0) 

Jo Jo Jo 

and {ui;)a{x') := (ui;)a(x', X 3 ) dx 3 . As a consequence, 

r^3 

{U;)aix) = {Ue)'^{x', 0 ) - da{U;) 3 {x', s) ds + (fi{x) 

Jo 

rx3 

da{Ue) 3 {x',.s)dsdX 3 - / daiUe) 3 ix', s) ds + ri^^{x), 

Jo 


(6.10) = {Ue)a{x') + 


1 rX 3 



0 JO 


where ||l?a||L2((a;\Ae)x(0,l)) ^ 0. 


26 


JEAN-FRANgOIS BABADJIAN, DUVAN HENAO 


On the other hand, for £^-a.e. x & Qj, let us define the sequences 

r^3 

9a{x\x^) := Xa;\Ae(a;') / 5a(Me)3(a:', s) ds, 

Jo 


9l{x') ■= Xa;\Ae(^') 


1 rxi 



9a(w£)3(a;',s)dsdx3. 


0 ^0 


From (16.9p and the a priori bound < M, we get 115^11^2(0^.) < C for some 

constant C > 0 independent of e. Therefore, up to a subsequence, ga weakly in 

L‘^{Qf) for some ga G T^(0/-). In addition, 5 ^ ^ ga weakly in where ga{x') := 

/o 5a(a;',a:3)dj:3^ _ 

Multiplying (16.101) by Xuj\Ae leads to 

{Ue)a{x)Xi^\A,ix') = iUe)aix')Xuj\AAx') + faix') - g"a{x) + faix), 
where || 5 alU 2 (Oj.) —>• 0. Passing to the limit as e —0 finally yields 
(1 - 0{x')){Ua{x) - Ua{x')) = ga{x') - ga{x), 
and according to the structure (|5.18l) of planar displacements, we deduce that 


- X 3 ^ (1 - 9{x))daUz{x) = ga{x) - ga{x). 


Since by assumption U 3 = 0 £^-a.e. in {0 < 1}, we get by locality of approximate gradients 
of SBV functions (see O Proposition 3.73 (c)]), that Vtt 3 = 0 £^-a.e. in {6 < 1}, hence 
ga{x) = ga{x'). As a consequence, XLj\Aeda{ue )3 = Dog% ^ ^sga = 0 weakly* in 
L2 (w; 77-1(0,1)). □ 

An alternative to the definition of as the orthogonal projection of Jue onto w x {0} 
is to consider its projection along certain almost-vertical oblique directions. Define the 
unit vectors 

^^ = ^(±1,0,1), 77± = ^(0,±1,1), 

and their rescaled versions 

Denote by 7 r^± (resp. 7r^±) : —)• the projection onto {xs = 0} parallel to the vector 

(resp. rjf), i.e., for x := (x',0) + (resp. x := (x',0) + trif), then 7r^±(x) := x' 
(resp. '/r^±(x) := x'). Finally, consider the set 

Ag . 7 r^+ (Aig (d (cUg X (—2,1))) U («Ate D (cug x (—2,1))) 

U TT^- (J„^ n (We; X (- 2 , 1))) U H (cJe X (- 2 , 1 ))) U P^. 

where uje ■= {x' G oj : dist(3:', Scu) > 2e}. Up to a subsequence, it can be assumed that 

XAe ^ 0 weakly* in L°°{oj) for some 0 G [0,1]). 






REDUCED MODELS FOR LINEARLY ELASTIC THIN FILMS 


27 


Using the decomposition 

= ^|e(i^uj2 + + ^\£{^ue)2 - (i^ue)3p + e^|(uuJll^ 

it is possible to prove that 

liminf f | (z^«e)3) | 

which shows that ||0||Li(tj) is controlled (up to a multiplicative constant) by the fracture 
energy in the bonding layer. The constant 1/8, however, is not optimal, since in order 
to obtain (16.61) that prefactor should not be present. In most situations [e.g. if the sets 
have uniformly bounded perimeters) it should be possible to obtain the optimal 
lower bound, but there are pathological cases (such as the microstructure Example 16.d|) 
where f 9 dx' is larger than the fracture energy on the left-hand side (because each vertical 
crack is counted twice in A^, which is defined as the union of all the oblique projections). 

Be it as it may, by including in the oblique projections of the cracks inside the 
bonding layer, one is able to obtain an optimal estimate for the elastic energy required by 
the body to accomodate the strain mismatch between the deformations in the film and in 
the rigid substrate. Before proving this final estimate, we need two preliminary technical 
results concerning sections of BH-functions along the oblique directions defined above. 
For £^-a.e. x' G oig and £^-a.e. t G (—2-\/2e, V^e), define the functions 

(ue)|±(t) := U£(^(x',0) :=Us(^{x',0)Tjf. 

Lemma 6.5. For C^-a.e. x' G Wg \ A^, we have 

(ugjfi G H^{—V2e,V2e) and G H^(—V2e,V2e), 

Se 'He 

with = 0, and 

S£ 'He 

X2, ^ (Ue)3(x',X3) G H^(0, 1). 


Proof. Let us denote by 

n^± := {C G m3 : C • = 0} 

the plane orthogonal to passing through the origin, and, for y G n^±, we define 

:= {t G M : ?/ + G H/}- 

According to slicing properties of functions of bounded deformations (see [H Theorem 
4.5]), we know that for "H^-a.e. y G n^±, the function 

1u^[y + t^f) ■ belongs to SBV‘^[Q^±), 

and its jump set is contained in 

{t G : y + G Jue}- 

Let us denote by N^± C n^± the exceptional set of zero measure on which the pre¬ 
vious properties fail. Since 7r^± are Lipschitz functions, it follows that the sets Z^± := 
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vr^±(A^^±) C oj are £^-negligible as well. Consequently, for all x' ^ uj^ \ Z^± (and thus for 
£ 2 .a.e. x' € uje)i we have that 

{ueT'± G SBV\-2V2e,V2e), 

se 

and its jump set is contained in 


{t G {-2\/2e,V2e) : (x',0) + G Jue}- 

By definition of the set Ag, if x' G \ Ag then (x', 0) + 0 Ju^ H [oos x (—2,1)] for all 

t G {—2^/2e, \/2e), and therefore (we)*! G H^{—2\^e,^/2e) for £^-a.e. x' G \ A^. In 

S£ 

addition since = 0 £^-a.e. in (— 2 \/ 2 e, —\/ 2 e), it follows that (—\/ 2 e) = 0 . 

se S£ 

The statement concerning the vectors can be proved in an analogous way. 


According again to slicing properties of functions of bounded deformations, we have 
that for C^-a.e. x' G w, the function X 3 i-A (ue) 3 (x',X 3 ) belongs to SBV‘^{Q,1), and its 
jump set is contained in {X 3 G (0,1) : {x',x^) G Jug}- a consequence, for C^-a.e. 
X £ io\ Ag, the function X 3 i-A {ue) 3 {x', X 3 ) belongs to H^{0, 1). □ 


The following technical result will be useful in the argument leading to a partial bulk 
energy lower bound. 

Lemma 6.6. Let G [0,27r), p := cos'd, q := sin??, a'nd defi'ne the unit vectors 

■■= ^{Tq,±pA)- 

For any matrix 'we have the decomposition 

|A|2 = |A^+• ^+p + |A^^-^~|2 + |Ar 7 +-77+12 +|Ar 7 ^-T 7 -| 2 -i(trA )2 

+P^«22 - ‘2^pqai2?‘ + + q^a22 + 2,pqai2\^ 

+2|(p2 - q^)ai2 + pq{a22 - aii)P + + (oii + 022 )^)- 

Proof. Let us define := ^+ A = {q,—p,0) so that {^+,^”,^ 0 } orthonormal 

basis of M^. Then the family 

{^+ ® ^+, r ® r, ^0 ® ^ 0 , o ^o), ^(r © ^o), © r)} 

defines an orthonormal basis of the set and Pythagoras Theorem ensures that 

I A|2 = IA : (^+ ® ^+)|2 + IA : (r © T )P + IA : (^q © ^o)l" 

+ 2| A : it © ^o)l" + 2| A : (r © ^o)l" + 2| A : {t © T)l" 

= \At • ^+1" + |Ar • r I" + |A^o • 

+ 2|A4+ • + 2|A^- • + 2|A^+ • 

The conclusion follows from a straightfoward computation of each term. □ 


We now prove a partial bulk energy lower bound. 
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Lemma 6.7. Assume that Xi, > fii,. Then (1 — 0)u^ = 0 C^-a.e. in u, and 


(6.11) liminf ilfo) > — liminf 

£^■0 2 £^>0 


j 

r^3 

Uaix)+ daiUe)3ix',s)ds 

/(a;\Ae)x(0,l) 

Jo 


dx. 


If in addition the sequences {da{u^)^)e>o o.re bounded in then 


hb 


liminf Ofe) >— / (1—0)|wrdx . 

£^•0 2 


Proof. Let us denote by 


£en{u^] 

) eei2('t^£) 

ei3(We) 

£ei2{Ue] 

) £e22{Ue) 

623 (We) 

eioiUe) 

623 ( We) 

e"^ 633 (u£) 

_h f 

2 Jn^ 

tr(A£)2 dx 

A Tb ^£ 

J rih 


the scaled strain so that 

According to Lemma lUlUl with the angle d = 0, we get that 
^6 hb 


(6.12) J,{u„nb)> 


'Ul, 


tr(A£)^ dx 


+ hb f + |Aer/+• T/+P + lA^r? ■ r] p] dx 

J Qi, 

> hb f [|-4.e^+ • ^+1^ + \AeC ■ + |Ae?7+ • r/+|^ + \Asr]~ ■ dx, 

JQf, 

since > fib. It remains to compute each of the four terms in the right hand side of the 
previous expression. Let us start with the first term. Changing variable x = (y',0) + 
(with dx = {y/2e)~^ dy'ds), and using Fubini’s Theorem, we get that 


• ^^l^dx > e- 


'fit, 


> e" 


L 


{LUe\Ae)x{-l,0) 
0 


•^j'Pdx 


[ / |Vit£((y',0)+ s^+)^+•^+pdsdy'. 

«/LJe\A£ J — y/ 2 € 


According to Lemma [631 since {u^y^ S H^{—V2s,V2s) and (wg)^^(—\/2e) = 0 for 
-a.e. y' G CO \ A^, we get that 


/ 

J 


2i£^’^-^^l^dx > e- 


I f 

JcjAA^ j — \/ 


e\Ae J-V2e 

rO 


ds 


[Ue{{y',0) + S^t) ■ ^t] 


ds dy' 


> e" 


[ / [ii£((2/',0)+ s^+)-^^Ids 

JuiAAr J-V2e hS 


\L 


0 Je\Ae 


dy' 


(^^£)i (y',0) + -(u£)3 (y',0) 


dy', 
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where ri“(-,0) denotes the lower trace of on a; x {0}. Using again Lemma 16.51 the 
function G 

4.5 (iv)], it follows that 


function (ue)^+ G iL^(—\/2e, \/2e) does not jump at t = 0. Thus according to [U Theorem 


+e =(«£)+ + e "H^-a.e. on w x { 0 }, 


and therefore, 


(6.13) 




'fit, 






(u,)+(y',0)--(u,)+(y',0) 


Analogously, we can show that 


(6.14) 

OjA.r-rpdx>ij 

f 

{Ue)tiy',0) - i(u£)(]'(y',0) 

(6.15) j 

dx > ^ ^ 

f 

t^e\Ae 

{Ue)t{y',0) + ^(U£);[(y',0) 

(6.16) J 

• r?“ pdx > ^ j 

f 

t^e\Ae 

(ue)^(y',0) - ^(u£);[(y',0) 


dy'. 


dy', 


dy', 


dy'. 


Summing up (16.131) . (|6.14l) . (|6.15l) . (16.161) and using (|6.12p leads to 
(6.17) Je{u^Mb)>^ ( (u£)+(y',0)(u£)+(y',0)dy'+ ^ / \{u^)^ {y', 0^ dy'. 

^ JluAAk ^ JujAAf- 


Since C Ag, Lemma 16.51 together with the fundamental Theorem of calculus yields, 

(6.18) [ \{Ue)3ix',X3)-{Ue)tix',0)\'^ dx<4: [ 1 633 (u^)dx < 

i(Lj\Ae)x(0,l) JUf 

In particular, (I6.17p . ()6.18p and the energy bound (16.51) ensure that 

f \{ueh?‘<^x <Ce^, 

i(a;e\Ae)x(0,l) 

which implies, letting e —)• 0, that (1 — 9)u3 = 0 C^-a.e. in uj. Therefore Lemma 16.41 shows 
that Xaj\Ae^a('^e)3 ^ 0 Weakly* in L^(a;; 77”^(0,1)). In addition, since C A^, we can 
use (16.91) and the fact that (ue)^ -A Ua strongly in L?‘{VLf), to obtain (|6.11l) . 

Assume now that the sequences {da{ue) 3 )e>o are bounded in Then the conver¬ 

gence of the planar gradient improves to XLj\Ae^aiue )3 0 weakly in L'^{Qf), and thus 
(| 6 . 11 l) gives 

liminf Je(w£,Ofe) > ^ / (1 
2 

since 6 is independent of X 3 . 


— 0)\ur dx = 


¥/» 


— 9)\u\‘^ dx' 


□ 
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